Spin relaxation in quantum dots with random spin-orbit coupling 
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We investigate the longitudinal spin relaxation arising due to spin-flip transitions accompanied by 
phonon emission in quantum dots where the strength of the Rashba spin-orbit coupling is a random 
function of the lateral (in-plane) coordinate on the spatial nanoscale. In this case the Rashba 
contribution to the spin-orbit coupling cannot be completely removed by applying a uniform external 
bias across the quantum dot plane. Due to the remnant random contribution, the spin relaxation 
rate cannot be decreased by more than two orders of magnitude even when the external bias fully 
compensates the regular part of the spin-orbit coupling. 



c 



C 

O 

o 



> 
o 
in 
o 
o 

in 
o 
-i— > 



o 
o 



X 



I. INTRODUCTION 

A quantum degree of freedom of an electron localized 
in a quantum dot (QD), i.e. its spin, is thought to be a 
useful tool for the realization of nanoscale devices that 
can be used for information processing. 1,2 The interac- 
tion of spins with the environment on the one hand allows 
necessary read and write procedures on the other hand, 
leads to losses of the information held by the system. 
Thus, a critical issue regarding the possibility to convert 
quantum dots into a hardware realization of quantum in- 
formation devices is the the ability to manipulate quickly 
the spins of electrons localized in quantum dots and to 
keep them in the desired states as long as necessary. The 
spin-orbit (SO) coupling in QDs plays a crucial role both 
for the spin manipulation and lifetime of the prepared 
spin states. For example, the SO coupling allows the ef- 
fective manipulation of spins by an external electric field 
due to the fact that the electric field influences the or- 
bital degrees of freedom, and, through the SO coupling, 
the spin states. The most interesting example of such a 
manipulation is the electric dipole spin resonance 3,4 , the 
effect that occurs when the electric field of the incident 
electromagnetic wave causes spin-flip transitions resonat- 
ing with the wave frequency. In this case the electric field 
is a much more efficient tool for manipulating the spins 
than the magnetic field. The spin states in quantum 
dots can be prepared and controlled by an external opti- 
cal field 5 too, thus allowing an optical realization of the 
write and read operations for applications in information 
technologies. 

The SO coupling in two-dimensional systems based 
on (OOl)-type structures is described by the sum of the 

Rashba 6,7 Hr = o:r (o~xky — a y k x ^j and Dresselhaus- 
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matrices, and k|| = — iV|| 



are the coupling constants, a are the Pauli 
(e/cfi)A|| is the in-plane mo- 
mentum of electron. Here e is the electron charge, and 
A is the vector-potential of the external field. and 
i?D terms arise due to the artificial macroscopic asym- 
metry of the structures and due to the microscopic inver- 



sion asymmetry of the unit cells, respectively. For holes 
the SO Hamiltonian is more sophisticated leading to a 
more complicated spectra of spin excitations and spin 
dynamics 9-13 . 

In GaAs/Al^Gai-^As structures 14,15 and Si-based 
transistors 16 , the SO coupling constants typically range 
from 1CP 10 to 1CP 9 eV-cm. It is important to men- 
tion that by applying an external bias across the quan- 
tum well, it is possible to manipulate the magnitude of 
a in InGaAs/InAlAs-based 17 and GaAs/AlAs-based 18-20 
systems and even change its sign by doping 21 . In the 
asymmetric Si/Sii-^Gea quantum wells investigated in 
Refs. 22,23 , where the Dresselhaus term is absent due to 
the unit cell inversion symmetry and the the band gap is 
relatively large, the doping-induced SO coupling is three 
orders of magnitude weaker than in zincblende systems 24 . 

The SO coupling not only provides an ability to ma- 
nipulate spins with an electric field, and thus, hope- 
fully, to design a spin transistor 25 , but also leads to spin 
relaxation. The Dyakonov-Perel' mechanism 26 , which 
requires random scattering of electrons by impurities 
and/or phonons, describes the spin relaxation of non- 
confined electrons in the bulk and in two-dimensional 
(2D) structures where the electron momentum A; is a rel- 
atively well defined quantum number. Here the orienta- 
tion of the spin precession axis in the SO field Hr + 
changes randomly through scattering events. The spin 
(7dp) an d momentum relaxation rates are inversely pro- 
portional to each other in this case. The resulting spin 



'k r, with r being the mo- 
rn 01 ) 



relaxation rate is 7 DP ' 

mentum relaxation time, and a depends on aJJ 1 " 1 ' and 
aR,. The spin relaxation can also arise from the random 
paths of electrons in regular systems of antidots 27 and 
the tunneling of holes in arrays of Si/Ge QDs 28 . Vari- 
ous mechanisms of SO coupling and spin relaxation were 
reviewed recently in Ref.29. 

For this reason it is necessary to understand the lim- 
its to the abilities to manipulate the electron spins in 
QDs, where, due to the localization, the electron mo- 
mentum is not a good quantum number, and thus the 
spin dynamics 30 , and hence the spin relaxation require a 
more sophisticated consideration. A possible " admixing" 
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mechanism of spin relaxation in QDs arises since SO cou- 
pling mixes the electron states with opposite spins, and, 
therefore, can cause spin-flip transitions by the emission 
a phonon 31 ~ 34 . This effect is a manifestation of an ef- 
fective spin-phonon coupling caused by the SO interac- 
tion. It has been shown, however, that other mechanisms 
can be important for the spin relaxation. For exam- 
ple, the interaction of the electron spin with the spins 
of nuclei in GaAs QDs can lead to spin relaxation even 
if no SO coupling in the Rashba or Dresselhaus form is 
present 35 . This interaction limits the maximal spin re- 
laxation time. In Si-based QDs the spin relaxation can 
arise due to the modulation of the electron g— factor by 
phonons 36 . However, theoretically considered 37-39 and 
recently observed 22,23 Rashba- and Dresselhaus-type SO 
coupling in 2D Si systems can lead to a more conventional 
"admixing" mechanism of the spin relaxation there. We 
note here that the randomness in SO coupling leads to 
a Gaussian rather than to an exponential decay of the 
spin polarization 40 as well as to a serious limitation of 
the operational modes of the proposed spin transistor 
devices 41 . 

Thus, the possibility to manipulate the spin states de- 
pends on the SO coupling strength, which, therefore, 
plays both a positive and negative role in the spin dy- 
namics of QDs. A solution to this dilemma can be found 
in the ability to change the strength of the SO coupling 
by switching it on and off when the desired spin states 
are produced and preserved, respectively. For example, 
the effects of the structural asymmetry of a quantum well 
can be strongly reduced 17,21 and/or the Rashba and the 
Dresselhaus terms can be effectively compensated by ap- 
plying an external bias 42 . If realized experimentally, the 
suggestion of Ref. 42 , is expected to open a route for a 
larger (up to 100 /im) size spintronic devices. 41 

Here we concentrate on the effect of random SO cou- 
pling on spin relaxation in QDs arising in systems where 
the artificial asymmetry and, in turn, the Rashba-type 
SO coupling is produced by one-sided doping. The ran- 
domness of the SO coupling in our model arises due to 
fluctuations in the dopant concentration. Because of the 
randomness, the SO coupling cannot be completely com- 
pensated by applying an external bias since the bias is 
uniform as a function of the in-plane coordinate. We 
show that even when the bias removes the regular part 
of the SO coupling, the spin relaxation rate is still fi- 
nite due to the residual random SO coupling. This effect 
limits a possible decrease of the spin relaxation rate to 
two orders of magnitude at most. As another example of 
the role of inhomogeneous SO coupling in QDs we note 
its influence on the localization effects 43 and on spin dy- 
namics of electrons in quantum rings interacting with the 
surrounding nuclei 44 . 

This paper is organized as follows. Section II presents 
a model of the random SO coupling in two-dimensional 
structures and provides us with the basis for the calcula- 
tions for QDs. In the third Section we present a theory of 
spin relaxation in QDs with random SO coupling, empha- 



sizing InGaAs-like structures. A summary of the results 
and a discussion of possible extensions of this work is 
given in the conclusions. 

II. UNIFORM AND RANDOM SO COUPLING 

In 2D systems where SO coupling arises due to one- 
sided doping, the qr parameter cannot be considered as 
a constant in space since fluctuations of the concentra- 
tion of the dopant cause its randomness. For doped Si 
and Ge bulk crystals the importance of randomness of 
the SO coupling was first understood by Mel'nikov and 
Rashba 45 . Below we consider the dopant fluctuations- 
induced Rashba field in lateral QDs, calculate the corre- 
sponding spin relaxation rate, and show that for applica- 
tions the randomness imposes important restrictions on 
the minimum spin- flip rate. 

A structure consisting of a 2D channel where the QDs 
are formed and a narrow dopant layer with a 2D con- 
centration n(r) of dopants with charge |e| is considered, 
with r = (r x ,r y ) being the 2D in-plane radius- vector. 
As the factor that determines the local strength of the 
SO coupling, we consider the z— component of the elec- 
tric field of the dopant ions E z (p) at a point with 2D 
coordinate p in the well symmetry plane. We assume 
that the SO coupling is a linear function of E z {p) with 
a R(p) = a so\e\E z (p), where aso is a system-dependent 
parameter that includes the influence of the electric field 
on the polarization of the electron wavefunction in a 
quantum well 48,49 . The z-component of the Coulomb 
field of the dopant ions is given by 

E Z { P ) = ^ J n(r)f(p,r)d 2 r 7 (1) 

where e is the dielectric constant. The integration is per- 
formed over the dopant layer. Function f(p, r) has the 
form: 

/(p ' r) ^ [(p-r) 2 ° + z 2 ] 3 / 2 - (2) 

Here p — {p x ,p y ) is two-dimensional vector characteriz- 
ing the electron positions in the conducting layer. We 
assume that the correlation function of the dopant con- 
centration is "white noise" in the form: 

((n(ri) -n)(n(r 2 ) -n)) =nS(r 1 - r 2 ). (3) 

Here n = (n(r)), and (...) stands for the average. 
The fluctuations are taken to obey Gaussian statis- 
tics, as is commonly considered in the theory of doped 
semiconductors 46 . With the increase of the distance 
between the layers, the fluctuations of ari(p) become 
smaller and smoother. The total asymmetry-induced 
Rashba parameter is a sum of a regular average (a) a and 
a random term comprising the zero mean p-independent 
contribution: 

«r(p) = (a) + a rnd (p), (4) 
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with (a)o — 2irasdne 2 / e. For the S(z — zq) — doping con- 
sidered here the magnitude of the random term is given 
^ V ( a md) = Q; Soe 2 V / ^/ e - z o 38 - By comparing this re- 
sult with (a)o one concludes that randomness of the SO 
coupling becomes important when y/ (a 2 nd ) ~ {a)o, that 



is at nz 2 ~ 0.1, which is a typical parameter value for 
zincblende quantum wells 50 ' 51 . We assume that by ap- 
plying an external bias V across the QD plane the mean 
asymmetry- induced term (a) a is decreasing as (a) a = 
(a)o (V c — V) /V c and that a critical bias V c , which is of 
the order of few Volts 17 ' 21 , fully compensates the mean 
SO coupling. 

The Hamiltonian of the coordinate-dependent SO cou- 
pling has to be written in the symmetrized Hermitian 
form: 



H 



R — 77 



<Jx |fcy,aR (p)} - Vy |fcx,Q!R(p)}] , (5) 

where jfcj,aR , (p)j stands for the anticommutator. 

Quantitatively, the spatial behavior of ckr(p) is char- 
acterized by the correlation function F aa (p 12 ) , where 
p 12 = (p l — p 2 ) /zo, of the random Rashba parameter 
determined by: 

(a R (p 1 )a R (p 2 ))) = (a 2 nd ) F aa (p 12 ) . (6) 

For the Gaussian fluctuations (see, for example 
Refs. 37,46 ), the correlation function in Eq.(6) necessary 
for describing spin relaxation of semiclassically moving 
itinerant electrons 37 is: 



Faa (P12) 



[r 2 + lf 2 f(r- Pl2 ) 2 + l 



3/2 



d 2 i 



(7) 



For the following investigations of the spatially ran- 
dom SO coupling in quantum systems 47 described by 
the Hamiltonian in Eq.(5) we need additional correlation 
functions of the random Rashba field and its derivatives, 
namely (i, j are the Cartesian indices x,y): 



«r(Pi) 



daR(p 2 ) 
OP? . 



"r(P2) 



\ d P f 
dan(Pi) da K (p 2 ) 



dp? d& 



1 2 \ Faa* (P12) 

("rnd) 1 , 

ZO 

I 2 \ F ai a (P12) 

("rnd/ " i 

z 

I 2 \ FctiOj (P12) fo s 
( a rnd) 72 • ( 8 ) 



The dimensionless correlators can be calculated, for ex- 
ample, as follows: 



F a a x (P12) 
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Other components of the correlation functions can be 
obtained from x < — > y permutations. We mention two 
important properties of the correlators concerning their 
behavior at small and large distances. The first property 
is that -Faa;— type correlators vanish at p l = p 2 . The 
second property is the fast decay of the long-distance 
asymptotic of the correlators given, for example, by: 



[r 2 + lf 2 (r-p 12 r + l 



3/2 



r y (r-pi 2 ) 2 



[r 2 + lf 2 (r-p 12 f + l 



5/2 



F aa (p 12 > 1) 



F aax (P12 » 1) 



~3 ' 
Pl2 

~4 : 
Pl2 



TTl I 1\ "12 "12 

f".", lPl2 » J-J ~ -75 ' 

Pl2 



(10) 



with the unit vector ni 2 = Pi 2 /pi 2 - The correlation func- 
tions shown in Fig.l decay at p 12 ~ 1, and therefore, es- 
tablish a general spatial nanoscale for the lateral and the 
z— axis directions. 

In zincblende-based structures the Dresselhaus SO 
term arising from the unit cell inversion asymmetry has 
to be added to the Rashba term. In the (001) structures 
the coupling parameter aj, 001 ' = a c (ir /w) 2 , where a c is 
the bulk Dresselhaus coupling constant (w 25 eVA 3 in 
GaAs and InAs [51]), and w is the quantum well width. 
We consider below a model where the random contribu- 
tion to the Rashba coupling (a 2 nd ) is much larger than 



, which is possible for sufficiently broad quan- 

J(001) 



turn wells, and, therefore, neglect the Hf^ ' term in the 
calculation procedure and subsequently discuss the role 
of this term in relation to our results. 

Another possibility is provided by (Oil) QWs, where 
the Dresselhaus term 8 has the form: 



H 



(on) 



a% 11] k y G z [l-(k 2 y -2k 



2 ^ W 2 TT 2 ] 



(11) 



d 2 r, (9) 



d 2 r. 



Here the z— axis is perpendicular to the QW plane and 
the in-plane axes are: x = [100] and y = [011], and 

a £> n ^ = a c{^/ w ) 2 /2. The Hp 11 *' coupling proportional 
to a z does not lead to the electron spin-flip, and, there- 
fore, the main contribution to the spin relaxation rate 
comes from the randomness. 

The strength of the Dresselhaus coupling can be con- 
siderably decreased by the strain induced in QWs by a 
mismatch in the parental lattices, which, therefore, can 
decrease the contribution of the Dresselhaus terms in the 
spin relaxation rate in (001) structures. 52 



III. LONGITUDINAL SPIN RELAXATION 
RATE 

First, we briefly review the influence of the random SO 
coupling on the spin relaxation of non-confined electrons 
in quantum wells. The spin precession rate and direction 
at a point p are determined by the random local cou- 
pling cur (p). As a result, the spin precession is random 
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even for a carrier moving straightforwardly and the total 
spin of the electron ensemble relaxes even if the regular 
term (a) a vanishes as a result of applying either an ex- 
ternal bias or equivalent doping at the other side of the 
conducting layer 37 . 

Now we consider a single-electron QD produced by ap- 
plying a lateral bias to a 2D electron gas (see Fig. 2) 
leading to electron localization. From the experimental 
point of view, we shall concentrate on Ino.5Gao.5As— like 
systems, where the ability to manipulate the SO coupling 
by a moderate bias applied across the structure has been 
clearly proven 17,21 . 

The spectrum of an electron in the QD, E n ^ z , is de- 
termined both by the confinement potential m*p 2 u> 2 /2 
(Fig. 2 shows two quantum dots on the random Rashba 
SO coupling background) and the uniform static mag- 
netic field B || z with the axial gauge A = [B,r]/2, 
respectively, leading to 



E n p z 



h 2 



: 2m*l 2 B 



+ 



(l + |4|+2n)±^B, 



2m*a 2 ' ' 2 

(12) 

where m* is the lateral electron effective mass, i z is the 
z— component of the orbital momentum, n is the radial 
quantum number, and the corresponding wave function 



spatial scale a = IqIb (Iq 



+ 1%) 1/4 , with l = ^h/2m*uj 
being the " oscillator" length in the absence of a magnetic 
field due to the lateral confinement forming the QD, hu> 
is the separation of the energy levels at B = 0, and \s) = 
I ]) 1 \{) is the spin state. With the increase of magnetic 
field, E n j z changes from the spectrum of a single-electron 
QD to the spectrum of a free electron in a magnetic field. 

We consider below the QD in only the orbital ground 
state with the function 



V>oo(p) 



1 



2ira 



exp(-p 2 /4a 2 ). 



(14) 



The spin- flip-transitions such as V'oo(p) IT) — * ^oo(p) U) 
without change of the orbital quantum numbers occur 
at a frequency lo = \g\ p B B /h and the transition en- 
ergy is released to an acoustic phonon with momentum 
Q = \d\ VbB/cxH, where c\ is the sound velocity for 
the given phonon branch A with A =|| and A =_L for 
the longitudinal and the transverse modes, respectively. 
The electron-phonon coupling Hamiltonian for acoustic 
phonons H c - P h has the form 



- P h 



-ph,P 



H 



-ph,D 



2ira 



eM-p 2 /4a 2 )L n ap/a)e u ^\s). 



(13) 

H ere ip is t he azimuthal angle, the magnetic length Is = 
\Jhc/ \e\ B, L„i z (p/a) is a Laguerre polynomial, and the 



arising due to the piezoeffect ^-H^-ph pj an< ^ due ^° ^ nc 
deformation potential ^-ff e -ph,£>^ • These terms have the 
form (assuming the crystal volume equal to one): 



H 



[A] 

e-ph,P 



q,e 



H. 



[A] 

e-ph,D 



q,e 



(15) 



where the summation is taken over the phonon mo- 
menta q, and polarization e. The propagation direction 
d = q/g, K14 is the strength of the piezointeraction, D 
is the deformation potential, £ is the crystal density, e 
is the dielectric constant, and fo^ qe is the phonon cre- 
ation operator. H^\ ph D is nonzero for the longitudinal 



phonon mode (de = 1) only. At finite temperature, the 
application of Fermi's golden rule yields the mean value 
of the spin relaxation rate = (T^ 1 ), where T\ is the 
relaxation time, with respect to the spin-flip as: 



<ri) 



2tt 



n B + 1] E (\9\ V B B) { (S'\ fl™ ph ^ gr (p) \s) 



(16) 



d.e 



where tib = [exp (\g\ p, B B/T) — 1] 1 is the Bose-Einstein 



phonon occupation factor. We consider below only the 
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zero temperature case for simplicity. The phonon den- 
sity of states is^ (e) = e 2 /2ir 2 (hc\) 3 , and (...) de stands 
for the average over the phonon directions and polariza- 
tions. The spin-flip matrix element of if _ p h is nonzero 
due to the admixing of the upper ip n £ z (p) \s) states to 
the ground state V> gr p) \s) , where \s) shows that spin is 
an approximate quantum number due to the SO cou- 
pling. Here we have neglected the spin splitting of the 
states g\i B B in comparison to the energy corresponding 
to the orbital degrees of freedom, hu> n i z . The approxi- 
mation is a reasonable one despite the large g— factor in 
Ino.5Gao.5As (g w 4.0) quantum wells, since a small ef- 
fective mass (to* = 0.04too) leads to gp B B /fou) n e z < 0.2. 
We mention here that due to the random coordinate de- 
pendence of the Hamiltonian H-r, there are no symmetry- 



related selection rules for n and l z in the matrix elements 
(s'| tp ni (p) ip m (p) \s) that are responsible for the 
formation of the spin relaxation rate, in contrast to the 
case of regular SO coupling 32 . Since the Hamiltonians 
P h.D nave drastically different phonon mo- 
mentum dependencies, the relative contribution of the 
deformation potential mechanism increases with an in- 
crease in the magnetic field 53 . Due to the large g- factor, 
this mechanism already dominates at B = 0.1 T. 

By using the approach suggested in Ref. 31 and tak- 
ing into account the Gaussian character of the fluctua- 
tions, after some algebraic transformations we obtain for 
the deformation potential contribution using the notation 
O=(0,0),l = (n l7 £f),2 = (n 2 ,^): 



( r i> - ( -^^f^ B5 J2^2 1 ~2 x ( 

n II 1,2 W m,<rna,<| 

'{K_{\)K*_{2)) (.9||(0; 1)^(0; 2)) ^ + {K-{1)K%{2)) (. 9 ||(0; l)$j[(2;0)) d + 
(^ + (l)^*(2))(. 9 ||(l;0)^(0;2)) d + (X + (l)^;(2))(.g||(l;0).^(2;0)) d . 

I 

The correlators are given by : with the matrix elements averaged over the disorder: 

(K_(l)iT(2)) = K vy + K xx -iK xy + iK yx , 

(K+{l)K* + {2)) = K yy + K xx + iK xy -iK yx , 

(K_(l)Kl(2)) = K yy -K xx -iK xy -iK yx , 

(K + {\)K*_{2)) = K yy -K xx + iK xy + iK yx , 



/" (2) 
^I'iP^n^iPliPj ^ (1) ^00(P 2 )^00(Pl)rf 2 Plrf 2 P2- (18) 

I 



The operators P^\P^ in Eq.(18) are defined according The matrix elements g\(n, £ z ; 0) of the phonon emis- 

to sion due to electron-phonon coupling are 

p(i) -_Ely -2i-^— p* 1 ) - Els. - 2i-^— (19) r n 

1% dp lx > ^ 1% Zt d Ply ' [ly) g x (n,£ z ;0) = J ^)C,,>)e^oo(p)d 3 r. (20) 

and act on the wavefunction V'oo(Pi) an d ^oo(P2)i re ~ where 4>o( z ) corresponds to the z— axis quantization and 

spectively. is assumed here for simplicity to be the rigid-walls wave 
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function tpo( z ) = \Z^/ W sin(7rz/u>). Since we consider 
broad quantum wells, where the size quantization energy 
is relatively small and the electron wavefunction weakly 
penetrates into the interfaces, the rigid- walls boundary 
conditions are justified. 

We now consider two applications of the model de- 
scribed above. First, we consider the effect of the fully 
compensated Rashba coupling (V = V c ) such that we 
are left only with the random contribution a rnc i(p) with 
(arnd(p)) = 0. Fig. 3 presents the results of numer- 
ical calculations (including a small contribution of the 
piezoeffect interaction) of the spin relaxation rate for 
Ino.5Gao.5As quantum dots as a function of the mag- 
netic field B using Eqs.(12)-(18) assuming D = 7.5 eV 
(Ref.[54]), K14 = -O.f Cm~ 2 (Ref.[55]), and sound veloc- 
ities C|| = 4.7 • 10 5 cm/s and c± = 2.8 • fO 5 cm/s for the 
longitudinal and transverse phonons, respectively (Ref. 
[55]). As can be seen, the spin-flip rate increases very 
rapidly with the applied field, analogous to the case of 
regular SO coupling 31 ' 53 . 

Second, we investigate the bias dependence of the spin 
relaxation rate (ri (V)) corresponding to the linear de- 
pendence of the mean Rashba parameter (a) a on the ap- 
plied bias. In this case, it follows from Eq.(18), that this 
dependence can be presented in the form: 

(ri(Q)) 
<ri (V)) 

™°» (21) 

(ri (v c )) + [(ri(o)) - (ri (v c ))\ ((v c - v) /v c ) 

Figure 4 shows the (ri(0)) / (T\ (V)) ratio for the two 
structures presented in Fig. 3. As can be seen in Fig. 4, 
the resulting (ri(V c )), (when there is full compensation 
((a) a = 0) ,) for z = 10 nm is approximately 60 times 
smaller than (ri(0)). The residual relaxation rate de- 
creases rapidly with the increase of the dopant-plane 
layer distance, corresponding to a decrease in the ran- 
dom fluctuations of the SO coupling. 

We note that the randomness of a rn d(p) leads to an in- 

/ 2\ 1/2 

homogeneous broadening or = ( (Ti — (ri)) ) of the 

spin-flip transitions making the relaxation rate depen- 
dent on the QD position. For the Gaussian fluctuations 
the broadening cr becomes of the order of the mean value 
of Ti at full compensation (V = V c ). The effect of in- 
homogeneous broadening can be seen clearly when the 
lateral confinement of the wave function a is consider- 
ably smaller than the spatial scale of the fluctuations of 
the random potential z . In this case each QD in the en- 
semble interacts mostly with the local Rashba field a(p ), 
where p is the position of the QD, that in this case varies 
weakly on the spatial scale of a. In the case of full com- 
pensation, the relaxation rate for each QD is proportional 
to the local Q^ nd (p ) leading to an inhomogeneous width 
of the transition of the order of (a 2 nd ) , being, therefore, 



of the order of (Ti) itself. Experimentally, this situation 
can be realized, for example, in strong magnetic fields 
B > 10 T, where the magnetic length Is < 10 nm. In 
the oCzo limit the ratio of the relaxation rates, which 
can be found from the discussion following Eq.(4), be- 
comes (T^Vc)) / (T, (0)) = (a 2 nd ) / (af = I/87rnz 2 . 

We now discuss the relative role of the random 
Rashba and Dresselhaus SO interactions since both ef- 
fects contribute to the relaxation rate at V — V c . In 
an Ino.5Gao.5As quantum dot of the width w =15 nm, 
as presented in Figs. 3 and 4, the Dresselhaus term 
a^ 01 ^ is of approximately 1.25 x 10~ 10 eVcm. At the 
same time, the one-side doping with n = 8 x 10 11 
cm~ 2 leads to (a) w 1CP 9 eVcm, giving in this case 
(ri(0)) / (T 1 (V c )) ~ 60. Therefore, the Dresselhaus and 
the random Rashba terms give similar contributions to 
the spin relaxation under these conditions. The ratio 
of these contributions depends on the applied magnetic 
field, and numerical analysis shows that the random part 
can dominate over the regular one. For (011) structures 
the random part gives the major contribution. 

IV. CONCLUSIONS 

In this analysis we have found that due to fluctuations 
of the SO coupling, the spin relaxation in QDs does not 
vanish even in the case of full compensation of the mean 
value of the SO coupling by an external bias. The resid- 
ual spin relaxation rate is of the order of few percent of 
the spin relaxation rate in the case of a non-compensated 
field and strongly depends on the system properties, the 
size of the dot, and the applied magnetic field. In the 
case of full compensation, the inhomogeneous line width 
of the spin-flip transition is of the order of the mean 
line frequency, and, therefore, one has a broad distribu- 
tion of the relaxation rates over the ensemble of QDs. 
The effect of the random Rashba SO coupling is compa- 
rable to or larger than the effect of the Dresselhaus cou- 
pling. The ratio of the contributions of these SO coupling 
mechanisms to the spin relaxation rate is compound- and 
structure-dependent, and shows that the importance of 
the randomness of the Rashba SO coupling. 

We concentrated here on the longitudinal spin relax- 
ation time Ti = T^ 1 that is related to the energy re- 
lease to the acoustic phonon bath in the spin-flip pro- 
cess. However, the spin ensemble dephasing time T27 
describing the motion of the spin component perpendic- 
ular to the magnetic field, and which may not be related 
to the energy transfer from the spin system to the lat- 
tice, is an important characteristic of a QD ensemble. 
As was recently found theoretically, the physics behind 
these two relaxation processes in systems of QDs cou- 
pled to phonons is similar 56 . Therefore, the randomness 
of the SO coupling considered above also influences the 
dephasing time T 2 . This is an interesting problem which 
will be investigated separately. The role of the magnetic 
field parallel to the plane 57 of the QD would be another 
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extension of the results presented here. 

In this paper we have discussed zincblende quantum 
dots only. However, some remarks concerning Si-based 
structures need to be given. If the admixing mechanism 
of spin relaxation in these structures becomes important, 
the arguments discussed above concerning the role of the 
randomness could be applied there to an even larger ex- 
tent due to the absence of the Dresselhaus term. A quan- 
titative analysis of this problem necessitates a compar- 
ison of the effectiveness of different mechanisms of spin 
relaxation in Si-based QDs, which is a separate interest- 



ing issue. 
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FIG. 1: Correlation functions of the random SO coupling parameter «r(p) and its derivatives: F aa (p 12 ) (solid line), F actx 
(p 12 ) (dashed line), F a:cay (p 12 ) (dash-dotted line); p 12 is parallel to the [110] axis. 
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200 nm 

FIG. 2: (Color online) Two quantum dots on a template with a random Rashba SO coupling. The arrows show schematically 
the spin orientations. The gates form the confining in-plane potential. The scale on the right represents the Rashba parameter 
in 10 -9 eVcm units, (a) = 0.5 x 10~ 9 eVcm, n — 2.5 x 10 11 cm -2 , and z = 10 nm. 
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FIG. 3: Logarithm of the spin relaxation rate as a function of applied magnetic field for different zq in the case of total 
compensation of the regular SO coupling by an external bias (V = V c ); zo = 10 nm (solid line), zo = 20 nm (dashed line), and 
lo — 15 nm. The relationship between the doping n and Rashba constant (a) corresponds to the experimental data of Ref. [21]. 




FIG. 4: The ratio of longitudinal spin relaxation rates (ri(0)} /(ri(V)) as a function of the applied bias that is gradually 
compensating the SO coupling. The applied magnetic field B = 1.5 T. 



